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We review some predictions of nonrelativistic effective field theories for heavy quark-antiquark systems with a 
small typical radius, r < 1/Aqcd- 



1. INTRODUCTION 

We are interested here in bound states com- 
posed only by heavy quarks (and gluons). In 
such cases, at least for the lowest states, the char- 
acteristic radius r (r being the qq distance) of 
the system is small, typically too small to probe 
the confinement effects, which arise at a scale 
1/Aqcd ^ r. In the following I will call such 
systems 'Coulombic' or ' quasi- Coulombic'. These 
situations are particularly interesting not only be- 
cause they have phenomenological relevance, but 
especially because they allow us to understand 
much more about QCD. In such cases both the 
mass scale m and the soft scale l/r are much 
bigger than Aqcd and thus still sit in the pertur- 
bative regime. Non-perturbative corrections exist 
but are not expected to be dominant. 

As it is apparent from the spectra, heavy 
quarkonia are non-relativistic systems. Thus, 
they may be described in first approximation us- 
ing a Schrodinger equation with a potential in- 
teraction. This amounts to saying that the heavy 
quark bound state is characterized by three en- 
ergy scales, hierarchically ordered by the quark 
velocity u 1: the quark mass m (hard scale), 
the momentum mv ~ l/r (soft scale (S)), and 
the binding energy mv"^ (ultrasoft scale (US)). In 
the Coulombic or quasi-Coulombic situation it is 
w ~ as. 

To address the multiscale dynamics of the 
heavy quark bound state, the concept of effec- 
tive field theory turns out to be not only helpful 



but actually necessary. QCD effective field theo- 
ries (EFT) with less and less degrees of freedom, 
can be obtained by systematically integrating out 
the scales above the energy we aim to describe. 
This procedure leads ultimately to a field theory 
derived quantum mechanical description of these 
systems. The corresponding EFT is called pN- 
RQCD [1]. Here, all the dynamical regimes are 
organized in a systematic expansion in the ratio 
of the different scales which is eventually an ex- 
pansion in V. 

In practice, by integrating out the hard scale 
NRQCD is obtained from QCD. After integrat- 
ing out the soft scale in NRQCD, pNRQCD is 
obtained. The Lagrangian of pNRQCD is or- 
ganized in powers of 1/m and r (multipole ex- 
pansion). The matching is done by comparing 
appropriate off-shell amplitudes in NRQCD and 
in pNRQCD, order by order in 1/to, as and or- 
der by order in the multipole expansion. The 
matching coefficients are non-analytic functions 
of r and have typically the following structure: 
V ~ V(r,p,Si,S2)(^lnmr + Bln^'r + C). 



2. pNRQCD for mv > Aqcd 

At the scale of the matching /z' {mv ^ n' ^ 
mz)^, Aqcd) we have still quarks and gluons. We 
denote by R = (xi -|-X2)/2 the center-of-mass 
of the qq system and by r = xi — X2 the rela- 
tive distance. The effective degrees of freedom 
are: qq states (that can be decomposed into a 
singlet S'(R, r, f) and an octet 0(R, r,t) under 
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Table 1 



Summary of the different kinematical situations. 



mv 




potential 


ultrasoft corrections 


> Aqcd 
» Aqcd 
» Aqcd 


> Aqcd 
~ Aqcd 
<. Aqcd 


perturbative 
perturbative 
perturbative + 
short-range nonpert. 


local condensates 
non-local condensates 
No US (if light quarks 
are not considered) 



o 



color transformations) with energy of the order 
of the next relevant scale, 0{h.QCD,'mv'^), and 
momentum^ p of order 0{mv), plus ultrasoft glu- 
ons (R, t) with energy and momentum of order 
0(Aqcd, Tii;^). Notice that all the gluon fields 
are multipole expanded. The Lagrangian is then 
an expansion in the small quantities p/m, 1/rm 
and in 0(AQCD,mt;^) x r. The pNRQCD La- 
grangian is given in this case at the leading order 
in the multipole expansion by[l]: 

£ = iv{s.(.*-^-r.M-i;i|;)s 

+.gVA(r)Tr{Otr -ES-hSV-EO} (1) 
+£'-^^^^Tr{0^r -EO-hOtOr -E} - ^Fj^^F""". 

All the gauge fields in Eq. (1) are evaluated in 
R and t. In particular E = E(R, t) and iDoO = 
idoO — glAoCR, t), O]. The quantities denoted by 
V are the matching coefficients. 

We call Vs and Vo the singlet and octet static 
matching potentials respectively. At the leading 
order in the multipole expansion, the singlet sec- 
tor of the Lagrangian gives rise to equations of 
motion of the Schrodinger type. The two last 
lines of (1) contain (apart from the Yang-Mills 
Lagrangian) retardation (or non-potential) effects 
that start at the NLO in the multipole expansion. 
At this order the non-potential effects come from 
the singlet-octet and octet-octet interactions me- 
diated by a ultrasoft chromoelectric field. Recall- 
ing that r ~ 1/mv and that the operators count 

^Although, for simplicity, we describe the matching be- 
tween NRQCD and pNRQCD as integrating out the soft 
scale, the relative momentum p of the quarks is still soft. 



like the next relevant scale, 0(mi;^, Aqcd), to the 
power of the dimension, it follows that each term 
in the pNRQCD Lagrangian has a definite power 
counting. This feature makes >CpNRQCD a suitable 
tool for bound state calculations: being interested 
in knowing the energy levels up to some power w", 
we just need to evaluate the contributions of this 
size from the Lagrangian. Prom the power count- 
ing e.g., it follows that the interaction of quarks 
with ultrasoft gluons (nonpotcntial or retardation 
effect) is suppressed in the Lagrangian by v with 
respect to the LO ( by gv if mv^ Aqcd)- 

In particular, pNRQCD provides us with the 
way of obtaining the matching potentials via the 
matching procedure at any order of the perturba- 
tive expansion[l]. In the EFT language the po- 
tential is defined upon the integration of all the 
scales up to the ultrasoft scale mv'^. From the 
matching to NRQCD in the situation Aqcd 
mv we can easily obtain the matching potential 
Vs at N^LL [1] 



Vsir) = -Cf 



av{r, IJ.') 



4tt 



7i5(4ai/?o + 2/3i) + ( y +47|)/32 



o:v{r, n) =as{r){l + (ai + 2jePo) 



16 7r2 



-1-02 



(2) 



where /3„ are the coefficients of the beta function 
(ofg is in the MS scheme), and ai and 02 were 
given in [2]. The Coulomb potential turns out to 
be sensitive to the ultrasoft scale but infrared fi- 
nite. The same happens with the other potentials 
(like Vo or the potentials that bear corrections 
in 1/m") that can equally be calculated via the 
matching procedure. The /x' scale dependence in 
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the potential is cancelled for any physical process 
by the contribution of the ultrasoft gluons that 
arc cutoff at the scale /i'. 

Then, there are two main situations [1], cf. 
Table 1. If mv mv^ > Aqcd, the sys- 
tem is described up to order a.^ by a poten- 
tial which is entirely accessible to perturba- 
tive QCD. Non-potential effects start at order 
Ofgln/i' [3,4]. We call Coulomhic this kind of 
system. Non-perturbative effects are of non- 
potential type and can be encoded into local (a 
la Voloshin Lcutwyler) or non-local condensates: 
they are suppressed by powers of Aqcd /mv^ and 
AQcn/mv respectively. If mv » Aqcd ^ mv^, 
the scale m,v can be still integrated out perturba- 
tively, giving rise to the Coulomb-type potential 
(2). Non-perturbative contributions to the po- 
tential arise when integrating out the scale Aqcd 
[1]. We call quasi- Coulomhic the systems where 
the non-perturbative piece of the potential can 
be considered small with respect to the Coulom- 
hic one and treated as a perturbation. Some lev- 
els of ti, the lowest level of bb may be considered 
Coulomhic systems, while the J/iIj, the rjc and 
the short-range hybrids may be considered quasi- 
Coulombic. The may be in a boundary situa- 
tion. As it is typical in an effective theory, only 
the actual calculation may confirm if the initial 
assumption about the physical system was appro- 
priate. 

For all these systems it is relevant to obtain 
a determination of the energy levels as accu- 
rate as possible in perturbation theory. In re- 
cent years our knowledge of heavy quarkonia 
in the framework of perturbative QCD has de- 
veloped significantly. On one hand computa- 
tions of new higher-order corrections to various 
physical quantities appeared e.g. [3,4,6]. On the 
other hand the discovery of the mechanism of 
the renormalon cancellation in the quarkonium 
spectrum led to a drastic improvement of the 
convergence of the perturbative expansion of the 
energy levels. As important applications up to 
date, the theory enabled precise determinations 
of the MS-mass of the bottom quark and of the 
top quark from (mainly) the energy spectra of 
the lowest-lying states. The main uncertainty 



comes, in the bottomonium case, from the (essen- 
tially) unknown non-perturbative contributions. 
These arc generally claimed to be around 100 
MeV and ultimately set the precision of the pre- 
diction. Within pNRQCD a complete and sys- 
tematic parametrization of perturbative and non- 
perturbative effects of the heavy quarkonium in- 
teraction can be performed (cf. Table 1 and 
[1]). The parametrization depends on the mutual 
relation between the scale of non-perturbative 
physics, Aqcd, and the other dynamical scales 
in the specific heavy quarkonium system under 
study. A way to determine it and the size and na- 
ture of the non-perturbative contributions, con- 
sists in establishing to which extent perturbative 
QCD can consistently and successfully describe 
the system. This is the problem we have ad- 
dressed in [5], which investigates the range of va- 
lidity of perturbative QCD on the heavy quarko- 
nium spectroscopy and consequently extracts up- 
per bounds to the non-perturbative contributions 
by comparing the perturbative predictions, at 
the current best accuracy, with the experimen- 
tal data. One of the main problems in having a 
consistent (i.e. convergent) perturbative expan- 
sion for the quarkonium energy levels has been 
for a long time its bad behaviour when expressed 
in terms of the pole mass. The reason has been 
understood in the renormalon language [7]: the 
pole mass and the QCD static potential, respec- 
tively, contain renormalon contributions of order 
Aqcd, which get cancelled in the total energy of 
a color singlet quark-antiquark system. The so- 
lution then consists in making explicit this can- 
cellation by substituting in the energy expansion 
a short-range mass (free from infrared ambigui- 
ties) for the pole mass. In order to make explicit 
the renormalon cancellation we will express the 
quarkonium energies in terms of the MS masses 
and rearrange the perturbative series in the so- 
called e expansion [8]. This will be the key in- 
gredient of our analysis. We note that according 
to a formal argument, the residual uncertainty of 
the perturbative expansion due to the next-to- 
leading renormalon contribution is estimated to 
be Aqcd x (ax Aqcd)^ for a bound-state X of size 
ax [1,5]. Provided this argument applies, in prin- 
ciple the predictions of perturbative QCD can be 
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Table 2 

The quantities Ej^ are the order contributions to the 
the experimental values of the J/ij) and T(15) masses, 
uncertainties in a^\Mz) only. All the other data refer to a^\Mz) = 0.1181 and to the msbar quark masses fixed on 
the central values. 



spectrum [5,8]. The c- and b-quark msbar masses are fixed on 
The uncertainties in the third and fourth columns refer to the 



State X 



Ex 



- Ex E' 



X 



E 



W 

X 



E 



X 



T(PS'i) 

T(l3Po) 
T(l3Pi) 

T(1-''P2) 
T(2'''S'i) 

T(2-''Po) 
T(2-''Pi) 

T(23P2) 

T(33S'i) 



9.460 

9.860 

9.893 

9.913 

10.023 

10.232 

10.255 

10.268 

10.355 

10.580 

6.4 ±0.4 



+75 



9.460 

9.995_62 
10.004t^^ 

+81 
-65 
-93 



10.012 
10.084; 



10.5481}^? 
10.564t^^^ 
10.578ti^^ 
10.645t?i« 



6.307^2 







-0.099] 
-0.061: 



-o.i35;^i 
-o.iii;^« 

-81 
^65 
-93 
h75 

-0.316;}«? 
-0.309;?°° 

-o.3io;?°3 

-0.290;?i8 
0.1±0.4 



0.866 
1.534 
1.564 
1.591 
1.618 
2.421 
2.472 
2.518 
2.472 
* 

0.675 



0.208 

0.101 

0.081 

0.063 

0.096 

-0.356 

-0.404 

-0.449 

-0.348 

* 

0.188 



0.006 

-0.021 

-0.022 

-0.022 

-0.010 

0.102 

0.116 

0.129 

0.140 

* 

0.017 



2.14 

1.08 

1.05 

1.034 

1.02 

0.778 

0.770 

0.762 

0.770 

* 

1.62 



0.286 
0.459 
0.468 
0.477 
0.486 
0.710 
0.726 
0.740 
0.726 
* 

0.334 



made precise down to this accuracy by sufSciently 
increasing the orders of the perturbative expan- 
sion. The complete analysis of the consistency 
between the whole level structure of perturbative 
QCD and the experimental data has been done in 
[5]. The major results of that analysis have been 
as follows (see Table 2). (1) Once the cancellation 
of the leading renormalon has been incorporated, 
the perturbative scries turns out to be convergent 
and to reproduce reasonably well the gross struc- 
ture of the bottomonium spectrum at least up to 
some of the n = 3 levels. (2) The constraints on 
non-perturbative and higher-order contributions 
to the bottomonium spectrum, set by the com- 
parison of the calculation with the experimental 
data, indicate that these are smaller than usually 
believed. Moreover, the perturbative predictions 
of the bottomonium spectrum up to n = 3 ap- 
pears to be in agreement with the experimental 
data (within the next-to-leading renormalon un- 
certainties and taken into account the error with 
which we know a^). 
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